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We calculate a tunneling time distribution by means of Nelson's quantum mechanics and investi- 
gate its statistical properties. The relationship between the average and deviation of tunneling time 
suggests the existence of "wave-particle duality" in the tunneling phenomena. 



I. INTRODUCTION 



CO 

o 

O 

CN 

F ^ ■ It was suggested that there is a time associated with the passage of a particle under a tunneling barrier, so-called 
| \ tunneling time 0. Actually, several authors have tried to measure the time experimentally However, there 

is no clear consensus about any definition of tunneling time admissible for everyone. There are several approaches 
and methods to estimate the tunneling time. For example, Wigner time 0, El 0] is based on the time evolution of 



. wave packet through the barrier, and the d elay time of the peak or the centroid is expressed by an energy derivative 
' of phase shift. Larmor time H, |j| 0, 0, 0] is obtained from the Larmor precession angle caused by a magnetic 
, field confined in the barrier region. The traversal time proposed by Biittiker and Landauer [lfll fLU llfij is defined 



t-H ■ by the analysis of transmission coefficient through a static barrier augmented by a small oscillation in the barrier 
height. This time is obtained by measuring or analyzing the effect of "clock" added on the tunneling barrier. The 
dwell time 0, fl6l Il7| is defined as the total probability of the particle within the barrier divided by the incident 

'J ' probability current. This is only applicable to the stationary state case. There are another type of methods based 

l„ T5„l : l ITol 1771 lorll tp ii. . — l l (oil lo^l lo^l 



on the motion of "particle paths" , for example, Bohmian mechanics [LsL ITU l2(i| , Feynman path integral pH 12^ . 125 
Nelson's quantum mechanics [271 123 . |29| . and so on. In these methods, the tunneling time is defined as the time spent 
by "particle paths" under the tunneling barrier. See Refs. |24l l25j and references therein for reviews of this problem. 

We think that Nelson's quantum mechanics |2d| have some characteristic properties to study the tunneling time 
as follows. Since this method is described by the real-time stochastic process, it enables us to describe the real-time 
evolution of individual events as the analogy of classical mechanics. We call such an event of evolution as a "sample 
path." Moreover, since a sample path has its own history, we can obtain information of time parameter, in particular, 
tT tunneling time. 

In this paper, first, we calculate a tunneling time distribution by means of Nelson's quantum mechanics |26| . 
The tunneling phenomena should occur quantum mechanically and accompany fluctuating properties. If a tunneling 
particle is described by a wave packet with distribution of wavenumbers of finite width, then one would have a 



distribution of tunneling time. The width of the distribution should be reduced by using the spatially wider wave 
packet or the sharper distribution of wavenumber wave packet. In such case there should still remain the fluctuation of 
tunneling time coming from quantum effect. In spite of such an argument, many approaches centered discussion only 
on the averaged value, because it is difficult to consider the tunneling time distribution by the use of the conventional 
frameworks of quantum mechanics. However, Nelson's quantum mechanics can afford to predict such distribution, 
because this method enables us to obtain an ensemble of various sample paths in tunneling barrier. Next, from the 
tunneling time distribution calculated by this method, we investigate the statistical properties of the distribution, 
such as the average and the deviation of tunneling time. We found that the relationship between the deviation 
and average suggests the existence of the "wave-particle duality" in the tunneling phenomena. Last, we discuss the 
"quantum-classical correspondence" by analyzing the "Planck constant dependence" of tunneling time distribution. 
In this paper, for simplicity, we analyze the tunneling phenomena with one-dimensional static rectangular potential 
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barrier with height Vq and width d as following, 



x < -d/2 (region I), 

V(x) = { V -d/2<x <d/2 (region II), 
x>d/2 (region III). 



(1) 



The initial Gaussian wave packet with variance Ax 2 , 



2n(Axy 



1/4 



exp 



4(Ax) 2 



■(p) (x- (x)) 



(2) 



is injected into the potential barrier, where (x), and (p) are an expectation value of position and momentum respec- 
tively. We use the natural unit m = % = 1, and perform numerical simulations of 100,000 sample paths. For the 
numerical simulation, we take [—1000, 1000] as the total space for the unit of 1/fco, and adopt these parameters as 
follows, (x) — — 500/fco, and (p) = ko- 



II. ESTIMATION OF THE TUNNELING TIME BASED ON THE NELSON'S QUANTUM MECHANICS 

In this section, we give a brief review of the Nelson's approach of quantum mechanics [2^| which plays a central role 
in this paper, and explain how to estimate the tunneling time. Nelson's quantum mechanics based on the real-time 
stochastic process, enables us to describe the quantum mechanics of a single particle in terminology of the "analog" 
of classical mechanics, i.e. the ensemble of sample paths. These sample paths are generated by the Ito type Langevin 
equation, 

dx(t) = [u(x(t), t) + v(x(t), t)]dt + dw(t), (3) 

where x(t) is a stochastic variable corresponding to the coordinate of the particle, and u(x(t),t) and v(x(t),t) are 
the osmotic velocity and the current velocity, respectively. The dw(t) is the Gaussian white noise with the statistical 
properties of 

(dw(t)) = 0, and (dw(t)dw(t)) = —dt, (4) 

m 

where (• • •) means the ensemble average with respect to the noise. In principle, the osmotic and the current velocities 
are given by solving coupled two equation, i.e. the kinetic equation and the "Nelson-Newton equation" . Nelson 
showed that, for the expectation value of the dynamical variable, e.g., x, p, the whole ensemble of sample paths gives 
us the same results as quantum mechanics in the ordinary approach |26| . Once the equivalence between Nelson's 
framework and the ordinary quantum mechanics is proved, it is convenient to use the relation 

u = Re — — hiip(x, t), and v = Im — — lnip(x,t), (5) 
m ox m Ox 

where ip is the solution of Schrodinger equation. Since an individual sample path has its own history, we obtain 
information on the time parameter, e.g. tunneling time definitely. 

Now using the Nelson's quantum mechanics, we estimate the tunneling time of a particle crossing over a potential 
barrier. First, we prepare an incident wave packet given by Eq.(j2Jl from the region I. Next, we solve the time-dependent 
Schrodinger equation. Last, using the relation Eq.(|SJ), we obtain the drift term of the Langevin equation ©, and 
calculate sample paths. Suppose a simulation of tunneling phenomena based on Eq.@, starting for t = — oo and 
ending t — oo. As we treat a wave packet satisfying the time-dependent Schrodinger equation, the wave packet is 
located region I initially and turns finally into two spatially separated wave packets which are regions I and III, 
respectively. Figure 1 shows a typical transmitted sample path calculated by Eq.©. Transmitted sample paths 
originate preferentially from the front of the initial wave packet as suggested by Imafuku et al.[27]]- Every transmitted 
sample path has its traversal time of barrier, i.e., the tunneling time which is described as 

*/ 

Q{xi(t))dt, {i=l,2,---,N), (6) 

o 

where Xi(t) is the i-th sample path, tr is the final time, the function <d(x) is unity for —d/2 < x < d/2 and zero 
otherwise. Collecting these events, we can construct a statistical distribution of tunneling time. 
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III. TUNNELING TIME DISTRIBUTION AND WAVE-PARTICLE DUALITY 



From the ensemble of sample paths, we define a distribution of tunneling time as follows, 

P(T)fr S *£2, (7) 

5n(r) is the number of sample paths with the tunneling time from t to r + St, and N is the total number of sample 
paths. 

Figures 3 and 4 show the average of tunneling time (r) and its deviation At versus potential barrier width for 
various width of wave packet Ax. We can see that in the case of Ax greater than 20/fco, behaviors of average and 
deviation are independent on Ax. However, in the case of Ax = 10/fco, the data deviate remarkably from the others in 
the region d > 10/fco- It is suggested that this effect originates from the wave packet spreading during the propagation. 

Indeed the spatial deviation of free wave packet at time t is J Ax 2 + 4rr X^ x i ■ Therefore, the wave packet spreads 

twofold at t = 2m(Ax) 2 /%. By taking account of the fixed time t — 500/fc 2 , in which the peak of initial wave packet 

arrives at the left edge of barrier, this occurs in the case of width Ax < ~ 16/fco. 

It seems that the numerical results of (r) are roughly similar to the WKB time twkb = md/hn, where k = 

\J 2m(Vo — Eq)/Ti 2 , expect for the case of Ax = 10/fco- However, let us examine much in detail those values in the 
thin region (d < 4/fco). Figure 5 is an enlarged copy of this part in Fig. 3. It has been shown that, in the opaque case, 
the numerical simulation gives almost same values of the WKB times [U . While we can see these features of 
opaque case in Figs. 3 and 5, the numerical values deviate from the WKB time in the translucent case characterised 
by small nd which is approximately less than 2. The WKB approximation is not proper in the latter case. This 
suggests the tunneling phenomena make a "phase transition" in a sense around nd ~ 2. 

Next we show in Fig. 6 the deviation At as a function of the average (r), which is calculated by changing the width 
d with fixed potential height Vq for several cases of incident energy Eq. For each case we can see a common feature 
characterized by the fact that At is proportional to (r) for nd < 2 and to a/ (t) for nd > 2. In order to check this 
feature quantitatively, we fit the tunneling time distribution using the Gamma distribution, 

P(r) = f3a+1T \ a+1) ^- T/0 («, ft r >0). (8) 

This distribution has the following statistical properties, 

(t) = P{a+ 1), and At 2 = f3 2 (a + 1). (9) 

If a is constant, the relation At oc (t) holds good, and this is a typical feature of such coherent phenomena as 
deviation versus average value of photon number in coherent photon state. On the other hand, if (3 is constant, the 
relation At oc yj (r) holds good, and this is a typical feature of such random phenomena as Poisson process and 
Brownian motion. Figure 7 is an example of fitting distribution. 

Figures 8 and 9 show the fitting parameters a and K 2 j3 as a function of nd, respectively. We see the nd dependence 
of fitting parameters is universal even for different potential heights. Note that k is kept a constant value for each 
case. Furthermore, we found that a is constant in the translucent region {nd < 2), and k 2 /3 is constant in the opaque 
region (nd > 2). Therefore, this assures the At-(t) relation as mentioned above. 

Anybody never doubts that the tunneling time is a pure quantum process. However, the tunneling time is not an 
observable in the quantum mechanics. Actually it can be closely connected to the time evolution of some specific 
observable of the tunneling system and it should be measured experimentally through the time dependence of the 
observable. Thus the statistical property of the tunneling time distribution should reflect the characteristic features 
of underlying quantum process. From the above discussion we think that the tunneling may occur coherently, or in 
mode of wave picture dominantly in the translucent case and randomly, or in mode of particle picture dominantly in 
the opaque case. The tunneling time distribution reveals the wave-particle duality in the tunneling phenomena. 



IV. QUANTUM-CLASSICAL CORRESPONDENCE 



We discuss the Planck constant dependence of tunneling time distribution. First, we introduce the parameter 
e (0 < e < 1) and the "Planck constant" fi = eh, and consider "Schrodinger equation" of a wave function tp(x,t), 



- d ~ 

tn—ip{x,t) 



h 2 d 2 

2m dx 2 



V(x) 



(10) 
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under the conditions with fixed values of m, Vb, and po(= Uko). This wave function describes a virtual quantum system 
with scaled Planck constant h. This equation is formally transformed into the ordinary Schrodinger equation 



ih^(X,T) = 



h 2 d 2 



W{X) 



$(X,T), (11) 



2m dX 2 

by a scale transformation as following: 

X = x/e, T = t/e, W{X) = V(eX), and T) = i>{eX, eT). (12) 

Here let us clear up the procedure of the simulation. We fix the parameters m, Vq, E = p 2 /2m, and d. We scale the 
space-time variables as X — x/e. T — t/e and the potential width as D = d/e. In this X — T reference frame de 
Broglie plane wave is scaled as e l ( kx ~ ut ) — > e i{kx~ujT) ^ w h ere — e fc an j q = 6UI . Then we perform the numerical 
simulation based on Eq. (|llf> and obtain the distribution of r in the X — T reference frame. Eventually we obtain the 
distribution of f , the tunneling time corresponding to a Planck constant %, by the scaling of f = er. 

Now we will examine how the average (f) = e(r) and the deviation Af = eAr depend on %. Let us start the 
simulation at a translucent case with e = 1, where the tunneling process proceeds in the wave mode. Then, we 
decrease the value of e gradually, the effective width of tunneling barrier becomes wider and the tunneling process 
should proceed in the particle mode, or in other words, quasi-classically. This procedure gives us the e-dependence of 
(f) shown in Fig. 10. We can see that (f) increases with a tendency to approach the WKB time as e is decreased in 
the region near e = 1, whereas it is almost independent of e in the region with smaller e. These tendencies may be 
understood as follows. As it is seen in Fig. 5, the simulated values of (r) fit almost to the WKB time, 



7WKB = ^2(V -£b) A (13) 



in the latter case. Thus, the tunneling time corresponding to fi is (f) = etwkb = y 2{v^-e ) ^" Actually note that the 

WKB time can be expressed only by classical quantities. On the other hand, in the translucent case, the simulated 
values are smaller than the WKB time and approach gradually it as the effective width of tunneling barrier becomes 
wider, that is, e is decreased. Moreover, these two tendencies cross each other at e ~ 1/2 which suggests the tunneling 
phenomena change their phase around hud ~ 2h. 

Next, let us consider the e-dependence of Af. In this simulation, we may guess that in the region near e = 1 (wave 
mode and At cx (t)) the deviation should be proportional to the average, 

Afoc(f), (14) 

whereas in the region with smaller e (particle mode and At cx yj (r)), it should behaves as 

Af cx yfey/Jf). (15) 

Moreover, the change from the former tendency to the latter one may occur around hud ~ 2h. 

We can see these tendencies in Fig.ll, which give another support of the idea of "wave-particle duality" in the 
tunneling phenomena. 



V. SUMMARY 



We calculate the tunneling time distribution by means of Nelson's quantum mechanics. From the resulting distri- 
bution, we derived the statistical properties of it, the average and deviation of tunneling time. First, we found that 
if an incident wave packet is so large as to look like plane-wave like, the Ax dependence of them is negligible. Next, 
fitting the data by Gamma distribution, we found that the shape of distribution is universally determined only by 
nd. Furthermore, by investigating the statistical properties of the distribution in two characteristic "translucent" and 
"opaque" regions roughly divided by nd ~ 2, we found that the "wave-particle duality" may be seen in the tunneling 
phenomena. Last, we consider the Planck constant dependence of tunneling time, introducing the parameter e as 
h = eh. Consequently, we found that the dependences of the average and the deviation suggest another support of 
the idea of "wave-particle duality" in tunneling phenomena. 

We are interested in the comparison between the tunneling time distribution based on the Nelson's quantum 
mechanics and that based on the other method, e.g., Bohmian mechanics. The detailed study of such a comparison 
is a subject in the near future. 
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FIG. 1: Typical transmitted sample path calculated by Eq.@. r shows the traversal time of barrier, i.e., the tunneling time. 
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FIG. 2: Numerical results of tunneling time distribution as a function of tunneling time with potential height Vo/Eq = 2.0. 
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FIG. 3: Numerical results of tunneling time average versus potential width with potential height of Vq/Eq = 2.0. 
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FIG. 4: Numerical results of tunneling time deviation versus potential width for potential height of Vq/Eq 
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FIG. 5: An enlarged copy of thin barrier region (d < 4/fco) in Fig. 3. 
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FIG. 6: Relationship between average and deviation of tunneling time with the width of wave packet Ax = 50/fco- Dotted line 
shows the relation At oc (t) . Dashed line shows the relation At oc U (r). 
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FIG. 7: Example of tunneling time distribution fitted by the Gamma distribution in the case of Vo/Eq = 2.0, d = 10/fco, <^x = 
50/fco- In this case, we choose the fitting parameters a = 11.3 and j3 = 0.79, using least-squares method. 




FIG. 8: Fitting parameter a versus nd for fixed k. 
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FIG. 9: Fitting parameter (3 multiplied by re 2 versus red for fixed re. 
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FIG. 10: e-dependence of tunneling time average with potential height Vq/Eo = 2.0 and width d = 1/fco (nd = 1). Dashed line 
shows the WKB time twkb = a/ 2(v "g "^y^- 
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FIG. 11: e-dependence of tunneling time deviation with potential height Vo/Eo = 2.0 and width d = 1/fco («d = 1). Dashed 
line shows the relation Af oc \[i\J (f) in the case of particle mode. 



